UNH-10-01 



Ground state energy of the interacting Bose gas 
in two dimensions: an explicit construction 

Silas R. Beane 

Department of Physics, University of New Hampshire, Durham, NH 03824-3568. 

Abstract 

The isotropic scattering phase shift is calculated for non-relativistic bosons interacting at low 
energies via an arbitrary finite-range potential in d spacetime dimensions. Scattering on a {d — 1)- 
dimensional torus is then considered, and the eigenvalue equation relating the energy levels on the 
torus to the scattering phase shift is derived. With this technology in hand, and focusing on the 
case of two spatial dimensions, a perturbative expansion is developed for the ground-state energy 
of N identical bosons which interact via an arbitrary finite-range potential in a finite area. The 
leading non- universal effects due to range corrections and three-body forces are included. It is then 
shown that the thermodynamic limit of the ground-state energy in a finite area can be taken in 
closed form to obtain the energy-per-particle in the low-density expansion, by explicitly summing 
the parts of the finite-area energy that diverge with powers of N. The leading and subleading 
finite-size corrections to the thermodynamic limit equation-of-state are also computed. Closed- 
form results -some well-known, others perhaps not- for two-dimensional lattice sums are included 
in an appendix. 

PACS numbers: 05.30.Jp,64.60.an,67.85.-d 
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I. INTRODUCTION 



The study of quantum mechanical scattering in a confined geometry is topical in several dis- 
tinct ways. Recently developed experimental techniques involving trapped ultracold atoms 
are able to alter spatial dimensionality , thus motivating an understanding of the quan- 
tum mechanical interactions among atoms as the number of spatial dimensions are contin- 
uously varied. Bose gases in two spatial dimensions are of particular interest as they are 
expected to have a complex phase structure which is quite distinct from their counterparts 
in three spatial dimensions p-18]. On the other hand, from the perspective of numerical 
simulation, scattering in a confined geometry is often a practical necessity. For instance, 
in lattice studies of quantum field theories, calculations are done in a four- dimensional Eu- 
clidean space time volume. For reasons of cost, the finite spatial and temporal extent of 
these volumes is currently not enormous as compared to the physical length scales that 
are characteristic of the particles and interactions that are simulated. Moreover, there are 
no-go theorems 0] for Euclidean quantum field theory that require a finite volume in order 
to extract information about hadronic interactions away from kinematic thresholds. The 
technology required to relate hadronic interactions to the finite-volume singularities that 
are measured on the lattice has been developed in Refs. 1CH15] and state-of-the-art Lattice 
QCD calculations have now measured the energy levels of up to twelve interacting pions and 



allowed a determination of the three-pion interaction 16]. Similarly, quantum Monte Carlo 
studies of many-body systems in nuclear and condensed matter physics are carried out in 
a finite volume or a finite area, and thus a detailed understanding of how the energy levels 
in the confined geometry map onto continuum physics is essential to controlled predictive 
power. 

The purpose of the present study is several- fold. First, we aim to present a general study 
of the ground-state energy of a system of N bosons interacting via the most general finite- 
range potential, confined to a finite area. This energy admits a perturbative expansion in 
the two-body coupling strength for the case of a weak repulsive interaction. As a necessary 
prelude to considering a confined geometry, we first review the subject of isotropic scatter- 
ing of identical bosons in d spacetime dimensions using effective field theory (EFT). It is 
assumed that the reader is aware of the advantages of EFT technology. We then present a 
general study of the relation between eigenenergies on a torus and continuum-limit isotropic 
scattering parameters, for any spacetime dimension. While the eigenvalue equation that 
we obtain is derived in the non-relativistic EFT, it is expected to be generally valid in an 
arbitrary quantum field theory up to corrections that are exponentially suppressed in the 
size of the geometric boundary. A general study along these lines in quantum field theory 
is quite involved and has been carried out only in four spacetime dimensions [ll[. The re- 
sults of that study demonstrated that boundary effects due to polarization are suppressed 
exponentially with spatial size and therefore the leading power law behavior can be found 
directly in the non-relativistic theory. Hence the leading effects are cap tured by the non- 



relativistic EFT, with relativistic effects appearing perturbatively [13J, [15]. In the case of two 



spatial dimensions, the exact two-body eigenvalue equation has been considered previously 
in the context of lattice QED simulations in three spacetime dimensions [17] . However, 
there is little discussion in the literature about the consequences of scale invariance in the 
confined geometry, and about the ground-state energy of the many-body system in a finite 
area. Moreover, to our knowledge, the closed-form results that exist for the relevant lattice 
sums in even spatial dimensions, which render this case a particularly interesting theoretical 
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playground, have not been noted previously. 

Following our derivation of the ground-state energy of a system of N bosons confined to 
a finite area, we demonstrate that the thermodynamic limit of this system may be taken 
explicitly, by summing the parts of the expansion that diverge in the large N limit. In the 
thermodynamic limit, the energy-per-particle admits a double perturbative expansion in the 
two-body coupling and the density. As a byproduct of taking the thermodynamic limit, we 
are able to compute finite-size corrections to the thermodynamic-limit formula. In addition, 
we trivially include the leading non-universal corrections due to three-body forces. Study 
of the weakly interacting Bose gas at zero temperature has a long history, beginning with 
the mean-field result of Schick 18], with subleading corrections computed in Refs [§[ 



There are some claims in the literature regarding discrepancies among the various studies. 
We will comment on these claims below. 

This paper is organized as follows. In section II we review low-energy non-relativistic 
scattering of bosons in the continuum. Using EFT we calculate the isotropic phase shift 
in an arbitrary number of spacetime dimensions. Section III considers low-energy non- 
relativistic scattering of bosons in a confined geometry, in particular onad- 1-dimensional 
torus. We obtain the exact eigenvalue equation which relates the energy levels on the torus to 
the two-body continuum-limit scattering parameters. In section IV we consider the ground- 
state energy of a system of Bosons confined to a finite area. We first develop perturbation 
theory on the d — 1-dimensional torus and recover the perturbative expansions of the two- 
body results found previously. We then focus on N bosons interacting via weak repulsive 
interactions in a finite area and give a general expression for the ground-state energy. In 
section V we demonstrate how to take the thermodynamic limit in order to recover the 
well-known low-density expansion, and we compare our results with other calculations. We 
also compute the leading and subleading finite-size corrections to the thermodynamic-limit 
energy density. Finally, in section VI we conclude. In two Appendices, we make use of 
some well-known exact results for even- dimensional lattice sums to derive some closed-form 
expressions that are useful for the case of two spatial dimensions, and we evaluate several 
sums involving the Catalan numbers, which are relevant for deriving the thermodynamic- 
limit equation-of-state. 

II. SCATTERING IN THE CONTINUUM 
A. Generalities 

Here we will review some basic EFT technology which will allow us to obtain a general 
expression for the isotropic scattering phase shift in any number of dimensions. If one is 
interested in low-energy scattering, an arbitrary interaction potential of finite range may be 
replaced by an infinite tower of contact operators, with coefficients to be determined either 
by matching to the full theory or by experiment. At low energies only a few of the contact 
operators will be important. The EFT of bosons 1 , destroyed by the field operator ip, which 
interact through contact interactions, has the following Lagrangian: 

C = V 1 fa + ^) V - %<AW - yV(^)V(^) - + • • • (!) 



1 For a review, see Ref. 



22\- 
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FIG. 1: Feynman diagrams that give the exact two-body scattering amplitude. The oval blob repre- 
sents the all-orders interaction derived from the Lagrangian. 



This Lagrangian, constrained by Galilean invariance, parity and time-reversal invariance, 
describes Bosons interacting at low-energies via an arbitrary finite-range potential. In prin- 
ciple, it is valid in any number of spacetime dimensions, d. The mass dimensions of the boson 
field and of the operator coefficients change with spacetime dimensions: i.e. [ip] = (d — l)/2, 
[C2n] — 2 — d — 2n and [D ] = 3 — 2d. While our focus in this paper is on d = 3, in our 
general discussion of two- and three-body interactions, we will keep d arbitrary as this will 
allow the reader to check our results against the well-known cases with d = 2 and d = 4. 
Throughout we use units with h = 1, however we will keep the boson mass, M, explicit. 

Consider 2 — > 2 scattering, with incoming momenta labeled pi, P2 and outgoing momenta 
labeled p'^p^. In the center-of-mass frame, p = Pi = — P2 , and the sum of Feynman 
diagrams, shown in fig. [TJ computed in the EFT gives the two-body scattering amplitude (22V 



24 



MP) = ~ ; E%C 2n > (2) 

where 

M fny f d D ~\ 1 

/oW = Y l2 ) J ( 27 r)^V-q 2 + ^ ' (3) 

and it is understood that the ultraviolet divergences in the EFT are regulated using di- 
mensional regularization (DR). In eq. (J3|) , /x and D are the DR scale and dimensionality, 
respectively, and e = d — D. A useful integral is: 



2 \2) J {2n) D - 1 ^ \p 2 -q 2 + i5 / 

- r5r-^T ( 3 -^-) (4) 
2 P [ P ' V 2 J (47t)^-D/2 ' l4J 



In what follows we will define the EFT coefficients in DR with MS. This choice is by no 



means generally appropriate [23l . |24|. However it is a convenient choice if no assumption is 
made about the relative size of the renormalized EFT coefficients. 

Now we should relate the scattering amplitude in the EFT, ^Qo), whose normalization is 
conventional and fixed to the Feynman diagram expansion, to the S-matrix. We will simply 
assume that the S-matrix element for isotropic (s-wave) scattering exists in an arbitrary 
number of spacetime dimensions. We then have generally 

e **to = 1 + iM{p)A 2 {p) , (5) 
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where Af(p) is a normalization factor that depends on d and is fixed by unitarity. Indeed 
combining eq. (J2J) and eq. (pjj) gives M{p) = — 2Im(Jo(p)) and one can parametrize the 
scattering amplitude by 



A 2 ( P ) 



Im(I (p)) [cot 8(p) - i] 



with 



cot S(p) 



Im(/ (p)) 



E C 2n P 



2/i 



- Re(J (p)) 



(6) 



(7) 



Bound states are present if there are poles on the positive imaginary momentum axis. That 
is if cot 5(27) = i with binding momentum 7 > 0. These expressions are valid for any d. 
In order to evaluate Io(p) it is convenient to consider even and odd spacetime dimensions 
separately. For d even the Gamma function has no poles and one finds 



hip) 



M 



m p 



d-3 



2(4vr)( d - 1 )/ 2 r (^fi) 



(8) 



As there is no divergence, the MS EFT coefficients do not run with \i in even spacetime 
dimensions. Hence the bare parameters are the renormalized parameters. For d odd, one 
finds 



Up) 



M 



P 



d-3 



2 (4 7 r)( d - 1 )/ 2 r(^fi) 



log - 



d- 1 



log 7T 



(9) 



where ipo{ n ) is the digamma function. Here there is a single logarithmic divergence, hidden 
in the 1/e pole. Hence in our scheme, at least one EFT coefficient runs with the scale \x. 
With these results in hand it is now straightforward to give the general expression for the 
isotropic phase shift in d spacetime dimensions: 



p d ~ 3 cot5(p) 



(4tt 



7rM 



-r 



d-l 



E^2n P 



2n 



P 



d-3 



2tt 



log 



where /7 is defined by equating the logarithm in eq. ( TTU1) with the content of the square 
brackets in eq. (Q. Note that this is an unrenormalized equation; the C2 n coefficients are 
bare parameters and there is a logarithmic divergence for odd spacetime dimensions. One 
must expand the right hand side of this equation for small momenta in order to renormalize 2 . 
It is noteworthy that the effective field theory seems not to exist for d > 3 and odd as the 
divergence is generated at leading order and yet requires a nominally suppressed operator 
for renormalization. 

The leading three-body diagram in the momentum expansion is shown in fig. [21 and the 
three-body scattering amplitude is given by 



^3 



- D 



(12) 



2 In the case of three spatial dimensions eq. (flU)) yields the familiar effective range expansion, 

pcotS(p) = — - + ^-r 3 p 2 + C(p 4 ) , 
a 3 2 



(11) 



with a 3 = MCo/(8tt) and r 3 = I67rC 2 /(MC^). 
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FIG. 2: Feynman diagram that gives the leading contribution to the three-body scattering amplitude. 



B. Two spatial dimensions 



In this section we consider the case d = 3 in some detail. This case is particularly interesting 
because of its analogy with renormalizable quantum field theories, and QCD in particular [25j 



261 ] . From our general formula, eq. f fTOj) . we find 

cot 6(p) = -log (4) T~\ + °iP 2 + °(P') ( 13 ) 



where 



MCM 8CM 

aM = ~ T - ; a2 = MClUT) ■ (14) 



Note that a 2 is a dimensionless coupling, and y\o^\ is the effective range. Neglecting 
range corrections, for ^(/z) of either sign, there is a bound state with binding momentum 
7 = /xexp(7r/2o!2(A t ))- In essence, this occurs because, regardless of the sign of the delta- 
function interaction, quantum effects generate an attractive logarithmic contribution to the 
effective potential which dominates at long distances. However, as we will see below, in the 
repulsive case, this pole is not physical. 

Many interesting properties in two spatial dimensions may be traced to scale invariance. 
Keeping only the leading EFT operator, the Hamiltonian may be written as 



H = rf 2 x 



(15) 



where the field and spatial coordinates have been rescaled by ip — > M 1 ^ 2 ^; x M~ x l 2 x. It 
is clear that classically there is no dimensionful parameter and indeed this Hamiltonian has 
a non-relativistic conformal invariance (Schrodinger invariance) (26[. This conformal invari- 
ance is broken logarithmically by quantum effects. Perhaps the most dramatic signature of 
this breaking of scale invariance is the vanishing of the scattering amplitude at zero energy, 
which follows from eqs. M and [£21 

The leading beta function in the EFT is 

A^CbM = ^C 2 M , (16) 

which may be integrated to give the familiar renormalization group evolution equation 

l-f£*2 {V) log f 
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It is clear from eq. (fTTj) that the attractive case, a 2 {lj) = ~\ a 2(^)\, corresponds to an 
asymptotically free coupling, while the repulsive case, a 2 {lj) = + 1 CK2 (m) I ? has a Landau pole 
and the coupling grows weaker in the infrared. We will focus largely on the latter case in 
what follows 3 . Note that the position of the "bound state" in the repulsive case coincides 
with the position of the Landau pole, which sets the cutoff scale of the EFT. This state is 
therefore unphysical. 

Below we will also make use of a more conventional 4 parametrization of the phase shift: 
coU(p) = -log (p 2 a 2 2 ) + cr 2 p 2 + 0(p 4 ) . (18) 

7T 

Here a 2 is the scattering length in two spatial dimensions. By matching with eq. [131 one 
finds a^ 1 = /iexp(7r/2a 2 (/i)), which in the repulsive case is the position of the Landau pole. 
Hence, in the repulsive case, a 2 x is the momentum cutoff scale. Therefore, from the point 
of view of the EFT, a 2 is a most unsuitable parameter for describing low-energy physics. Of 
course, while the parameter a 2 is expected to be very small as compared to physical scales, 
its effect is enhanced as it appears in the argument of the logarithm. 



III. SCATTERING IN A CONFINED GEOMETRY 

A. Eigenvalue equation 

With the scattering theory that we have developed we may now find the eigenvalue equation 
in a confined geometry with periodic boundary conditions. Specifically, we will consider 
scattering on what is topologically the (d — l)-dimensional torus, X^ -1 = SjL x <SL x • • • x 
<S(V-iv m the confined geometry, all bound and scattering states appear as poles of the S- 
matrix, or scattering amplitude, A 2 (p). Hence, from eq. fl2j) we have the eigenvalue equation 
M{p)~ l = 0, or 

scb^ = ^w . tfw = fz^E^. da) 

where we have chosen to define the sum with a sharp cutoff (d — 2 is ultraviolet finite). The 
sum is over k = 2irn/L where n G Z d_1 = (rii,n 2 , . . . ,rid-i) takes all integer values. It is 
therefore convenient to write 

# GO = ^- 2 L 3 ~ d £ -^— 2 , (20) 

where q = pL/2n and therefore A = 27rA n /L. As the EFT coefficients are defined in DR, 
we can write the eigenvalue equation as 

1 



MI { o DR \p)) = IUP) ~ Re(/ {A} (p)) . (21) 



E C 2n P 



,2ri 



For a recent discussion of the implications of scale invariance for many-boson systems in the case of an 



attractive coupling, see Ref. 27 1. 



4 With 0,2 = ae 7 /2 and 02 = a 2 /27r, this parametrization coincides with a hard-disk potential of radius 
a [ljj]. As we will discuss below, there appears to be some confusion in the literature as regards the 
distinction between 0,2 and a. 
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Here we have subtracted off the real part of the loop integral using different schemes on 
the two sides of the equation; the integral on the left is evaluated using DR and the one on 
the right is evaluated with a sharp cutoff A. The purpose of this procedure is to leave the 
renormalization of the EFT coefficients, which is of course an ultraviolet effect, unchanged 
while defining the integer sums using an integer cutoff. We then have via eq. (j7|) our general 
form for the eigenvalue equation 



cot S(p) 



It is straightforward to find 



1 



Im(io(p)) 



(?) - Re(/ {A} (p)) 



4 a} (p) 



M 



A 



d-l 



(4vr)^r(^ 1 ) (d-l)p 2 



iTx 1, 



d-l d+1 A : 



2 ' p 2 



where 2^1 is the hypergeometric function. 

The exact eigenvalue equation in d spacetime dimensions can be written as 



(22) 



(23) 



q d ~ 3 cot 5(p)=T 



d-l 



7T 



d+1 
2 



n" 



+ 



2A 



d-l 



ir(d — l)q- 



-Re 



iH 1 



d — 1 d 



1 A 

1 '• 

Q 



(24) 



where it is understood that A n — > 00 on the right hand side. This equation gives the 
location of all of the energy-eigenstates on the (d — l)-dimensional torus, including the 
bound states (with p 2 < 0). The binding momentum in the confined geometry reduces to 
7 as L 00. While the derivation given above is valid within the radius of convergence of 
the non-relativistic EFT, this eigenvalue equation is expected to be valid for an arbitrary 
quantum field theory in d dimensions up to corrections that are exponentially suppressed in 
the boundary size, L. One readily checks that eq. EH gives the familiar eigenvalue equations 
for d = 2 [28| and d = 4 [l0-12]| and is in agreement with Ref. [l7[ for d = 3. 



B. Two spatial dimensions 



In a finite area, the energy levels for the two-particle system follow from the eigenvalue 
equation, eq. 



cot S(p) 



7T 



*'§)+**(i) 



where 



5 2 W = J2~ j - 2vrlogA n . 



(25) 



(26) 



Using the results derived in Appendix II, this integer sum can be expressed as 

,2 



l 00 t-iY 

SAD = - ? +v> -^gi-Lfcii 



7f 



(2£ + l) 



(27) 
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where ?/>o is the digamma function, and V2 is defined below. 

We can now combine our low-energy expansion, eq. (Tl3|) . with the eigenvalue equation, 
eq. (P5D, to find 

Using the renormalization group equation, eq. [T71 we then have 

cot5'(p) = f^) , (29) 

where 

cot5'(p) = -— + a 2 p 2 + 0(p 4 ) . (30) 

and «2 = «2(2vr /L). We see that in the eigenvalue equation, the logarithms of the energy 
cancel, and the scale of the coupling is fixed to 27r/L, the most infrared scale in the EFT 5 . 
Therefore as one approaches the continuum limit, the repulsive theory is at weak coupling 
and the attractive theory is at strong coupling. 



C. Weak coupling expansion 



When the two-body interaction is repulsive, the eigenvalue equation, eq. [301 allows a weak 
coupling expansion of the energy eigenvalues in the coupling a 2 . For the purpose of obtaining 
this expansion, it is convenient to rewrite the eigenvalue equation in terms of the scale- 
invariant momentum q = pL/(27r). If one expands the energy in terms of the coupling one 



can write q 2 = q 2 . + tq\ + e 2 q 2 + . . ., and the eigenvalue equation becomes 



1 

« 2 



+ 



02 (2tt) 
L 2 



+ eql + ...) + ... 



rS 2 (?) 



(31) 



Note that in this expression, the range corrections break the scale invariance with power 
law dependence on L. Indeed, in the presence of the range corrections, one has a double 
expansion in a 2 and in 1/L 2 . It is now straightforward to compute the energy perturbatively 
by expanding eq. [31] in powers of e and matching. 
With q = (0, 0) one finds the ground-state energy 



E n 



Aa 2 
ML 2 



+ 0{a\ 



+ 



16 «2 CJ 2 
ML 4 



1 + 0(a 2 )) + 0{L- & ) 



(32) 



5 The prime on the phase shift indicates that the part of the scattering amplitude that is logarithmic in 
energy is removed. This is a consequence of the confined geometry. 
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where we have included the leading range corrections and 

An 

V 2 = - 27rl °S A « = 47rlog(eivr-3r(f)) = 2.5850; 

n^O 11 

°° 1 27T 2 °° 1 7T 3 

7,4 = J2 —± = — c = 6 - 0268 ' p 6 = Yl ^ = y c{3) = 4-6589 ' (33) 

n^O n^O 

where 7 is Euler's constant, C is Catalan's constant and ((n) is the Riemann zeta function 6 . 
Note that one can use the renormalization group to sum the terms of the form 0(a 2 n V 2n ). 
One finds, for instance, for the universal part, 



4aS 



ML 2 



1 \ 2 / ./ x 3 

2 

7T 2 I \ 7T Z 



(34) 



where = a2(27i/Lg) with ^ = e 7 ' 2 ^ 27r ' ) = 1.5089. In this expression, the corrections to 
leading order begin at (9(a' 2 3 ). This freedom to change the scale at which the coupling 
constant is evaluated will be essential when we consider the many-body problem below. 

IV. N BOSON ENERGY LEVELS IN A FINITE AREA 

A. Perturbation theory for two identical bosons 

Exact eigenvalue equations for energy levels of N bosons (with N > 2) in a confined geometry 
are not available in the EFT of contact operators in closed analytic form. Hence it is worth 
asking whether the energy eigenvalues of the N-body problem admit sensible perturbative 
expansions about the free particle energy. It is straightforward to approach this problem 
using time-independent (Rayleigh-Schrodinger) perturbation theory. We will first consider 
the two-body case. Consider the two-body coordinate-space potential, 

V(r u r 2 ) = V 2S d - 1 (ri-r 2 ) , (35) 

where rj 2 is the two-body pseudopotential, an energy- dependent function, which is deter- 
mined by requiring that the potential given by eq. [33 reproduce the two-body scattering 
amplitude, eq. El The single-particle eigenfunctions in the confined geometry are 

<r|P> = JJ^rn^- (36) 
The momentum-space two-body potential in the center-of-mass system is then, 

V p , p > = (-p,p|V|-p',p'> = , (37) 

where |— p, p) are the two-body unperturbed eigenstates with energy E® = p 2 /M = 
(27rn/L) 2 /M. The perturbative expansion of the energy for momentum level n is given 



6 These results are derived in Appendix I. 
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by: 



47r 2 n 2 772 

-C/n = -^^r + 



ML 2 L d ~ x 



1 + 



r/ 2 M 



(2vr) 2 L rf - 3 



+ 



r/ 2 M 



(2n) 2 L d - 



E 



n 2 — m 2 



E 
E ■ 



n 2 — m 2 



n 2 — m 2 ) 2 







^2 
L d-3 



.(38) 



Hence for three spatial dimensions we have the nice perturbative sequence 
{l/L 2 , l/L 3 , l/L 4 , . . .}. However, in two spatial dimensions we have {l/L 2 , l/L 2 , l/L 2 , . . .} 
with an energy independent two-body pseudopotential, and therefore there is no perturba- 
tive expansion in l/L about the free energy, as expected on the basis of scale invariance. 
However, there is, of course, an expansion in r} 2 itself. 

It is now straightforward to recover the perturbative expansion of the two-body ground 
state energy in the case of two spatial dimensions. Here one finds 



V2 



4:a 2 



(l + \a 2 a 2 (r + r 2 )) 



(39) 



where the momenta have been written as arising from a Hermitian operator. In the relation 
between the pseudopotential and the amplitude, the minus sign is from moving from the 
Lagrangian to the Hamiltonian and 1/N! is a combinatoric factor for N identical bosons. In 
order to deal with the divergent sums in eq. [381 we renormalize as in the exact case (eq. I2~TI) . 
and replace, for instance, the leading divergent sum with 



M 



2(2vr) 2 L rf - 



A„ 

E 



- Re(J {A} (p)) + Re(I^ (p)) 



{DR} 



111- 



(40) 



With d = 3 and n = 0, this expression becomes 



M 



2(2tt) 



(41) 



The scheme dependent part of the DR integral then defines the running coupling a 2 (/i). 
Hence, inserting eq. [39] in eq. [381 an d noting that the additional /i-dependent piece in 
eq. ST] sets the scale of the coupling a 2 to 2n/L as in the exact case considered above, one 
immediately recovers eq. [321 including the leading range corrections. We emphasize that 
the language of pseudopotentials used here provides convenient bookkeeping in developing 
perturbation theory, however it is not essential. 



B. Perturbation theory for N identical bosons 



In this section, we generalize the perturbative expansion of the ground-state energy to a 
system with N identical bosons. The coordinate-space potential for the N-body system is 

N N 

K( ri ,...,r N ) = m ^^(Ti-Tj) + % S d - 1 (T i -T k )6 d - 1 (T j -T k ) + ... , (42) 

i<j i<j<k 
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where the dots denote higher-body operators. We have 



Vs = ~A 3 = l -D, , (43) 

where we have used eq. [T2j It is straightforward but unpleasant to find the ground-state 
energy of the N boson system in perturbation theory with this potential. In the case of three 
spatial dimensions, this has been worked out up to order l/L 7 [l3-15||. The calculation in 



two spatial dimensions is essentially identical, as the combinatoric factors for the ground- 
state level are independent of spatial dimension, and therefore the dependence on spatial 
dimensionality resides entirely in the coupling constant and the geometric constants. 



C. The ground-state energy 



In the case of two spatial dimensions one finds the ground-state energy 



En 



Aa 2 
ML 2 



a 2 



7T 



v 2 + 



Ol 2 



7T 



+ 



02\ 

-6 

+ 



K 3 + 3 



N 



V 2 V A 



2 ) 

N 

2 



V 2 — 

1 2 



- 12 



P 6 -24 



16 a\ a 2 fN 
ML 4 



W 2 V A + 51P 6 - 2 

1 A)/N 



P fi - 300 



P fi -90 



V 2 V i + 2Q Q + n -v 6 



v 4 

N 

2 



o(<4) 



+ 



2 J L 4 6 

where (™)=n!/(n - k)\/k\, the V 2s are available in eq. [331 an d 

1 



(44) 



Qo 



He, 



EE 

n^O m^O 

A„ 

EE 

m^O n 



n 2 m 2 (n 2 + m 2 + (n + m) 



16.3059 



- nV 4 log A r 



-1.8414 



rrr 1 n I m- 



n + m 



(45) 



(46) 



These double lattice sums have been evaluated numerically. This expression for the ground- 
state energy is complete to 0(o4), and includes the leading non-universal effects due to 
range corrections and three-body forces. Expanding out the binomial coefficients gives, for 
the universal piece, 



E 



Aa 2 /N 



ML 2 \ 2 



«2 

V 2 



V2 + O) 2 ^ 2 + ( 2N - 5 )^) 



(5) ( V2 + ( 2N - 7 )^4 + (5N 2 - 41N + 63)P 6 + 8(N-2)(2Q + fto)) 
+ (p) 4 (^2 4 - 6P 2 2 P 4 + (4 + N-N 2 )P 2 + 4(27- 15N + N 2 )P 2 P 6 



+ (14N 3 - 227N 2 + 919N - 1043)P 8 + . . . I + <D(a\ 



(47) 
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Here we have included a part of the O^a^) contribution [15| for reasons that will become 
clear in the next section. The dots represent double and triple lattice sums that appear at 
(HI, [29[ and which we do not consider here. 
As «2 = ct2(27r/L), it is clear that this expansion is valid only for repulsive coupling, 
which is small in the infrared. The expansion is expected to be valid for (a^N <C 1 7 . The 
chemical potential and pressure are readily available from the ground-state energy via the 
formulas 



dEn 



dN 



_ 1 dE 

•p _ 



2L dL 



N 



(48) 



By inspection of the binomial coefficients in eq. HUone sees that the leading effects from 
three-body forces enter at 0(a\) and 0(0.2) through two-body interactions. Other three- 
body effects enter by way of the contact operator in the Lagrangian and appear at the same 
order as effective range corrections: that is, they are suppressed by L~ 2 as compared to the 
leading two-body contributions, treated to all orders. This is of course a consequence of scale 



invariance. It is worth noting the contrast with the case of three spatial dimensions [13 



15] . There, the two-body contributions to the three-body force are logarithmically divergent 
in the ultraviolet and are renormalized by the three-body force contact operator in the 
Lagrangian. Both effects appear at 0(L~~ 6 ) in the expansion of the ground-state energy. 
In two spatial dimensions, scale invariance ensures that the sums in eqs. H5] and H6] are 
convergent as there is no scale-invariant counterterm available. Moreover, this ultraviolet 
finiteness persists to arbitrary order in 0:2- 



V. THE THERMODYNAMIC LIMIT AND THE DENSITY EXPANSION 

A. The Lee-Huang- Yang strategy 

The underlying scale invariance of the two-dimensional system allows a great deal to be 
learned about the thermodynamic limit directly from from the expression for the N-body 
ground-state energy in a finite area. By thermodynamic limit we intend the limit where N 
and L are taken to infinity with the density, p = N/L 2 , held fixed. Our strategy will be as 
follows. First we will use the renormalization group equation for the coupling ct2 to change 
the scale at which the coupling is evaluated to a quantity that is finite in the thermodynamic 
limit. We will then rearrange the expansion of the energy into series determined according 



to the degree of divergence with large N 30]. These series must sum to finite quantities in 



the thermodynamic limit. All quantities that are finite in this limit are kept. We will see 
that this strategy will enable us to constrain the form of the low-density expansion of the 
energy density of the two-dimensional Bose gas. Moreover, we will see that, unlike in the 
case of three spatial dimensions, the series that are most divergent with N can be explicitly 
evaluated. 



7 For a more accurate measure of the regime of applicability of the expansion, see Appendix II. 
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B. Universality and broken scale invariance 



As the coupling a 2 in eq. H7]is evaluated at the far infrared scale 27r/L, a change of scale 
is required before performing the thermodynamic limit. Consider a change of scale to v — 
27r\Ap, where A is a number which represents the inherent ambiguity in the choice of scale. 
With this choice, a 2 is finite in the thermodynamic limit, and constitutes a small parameter 
in the low-density limit (assuming that A is a number of order unity.) Using eq. [T71 we can 
then reexpress the energy as 



En 



Aa' /N 



ML 2 \ 2 

„/ \ 2 



+ 



V 2 + vrlog(NA) 
V\ + (2N-5)P 4 + Trlog(NA) (w 2 + Trlog(NA))) 



- {^\ (V\ + (2N-7)P 2 P 4 + (5N 2 - 41N + 63)P 6 + 8(N-2)(2Q + K 
+ 7T log (NA) (3 (Vl + (2N - 5) Va) + vr log (NA) (3V 2 + tt log (NA) 



+ l^f ) (^2 - + (4 + N-N 2 )P 2 + 4(27 - 15N + N 2 )P 2 P 6 



+ (14N 3 - 227N 2 + 919N - 1043)P 8 + . . . ) + 0( 



a 



(49) 



where now a 2 = a 2 {v). This expression is independent of A up to 0(a' 2 ) corrections. The 
strategy is to rearrange the expansion according to the maximum powers of N that appear 
at each order in a' 2 . We can then re- write eq. H9]as the energy-per-particle: 



Eq 

N 7 



2a^ 
M 



1 

1? 



N 



G 



1 

N2 



71 log (N\)U + 1 



^) (V 2 + vrlog(NA) 



+ 



a-. 



TT" 



tt log (NA) ( 2V 2 + vrlog(NA) 



0(a> 2 3 ) 



where 



G{z) 
H{z) 

X(z) 



2z 2 V 4 - 5z 3 V 6 + Uz 4 V 8 + 0(z 5 ) 
-Qz 3 V 4 + 20z 4 P 6 - 70z 5 V 8 + 0(z 6 ) 

-z 3 ( 2V 2 V A - 41P 6 + 8(2Q + TZ ) ) 

+ z 4 (4V 4 V 6 - V\ + 227P 8 + . . .) + 0(z 5 ) 



(50) 



(51) 
(52) 

(53) 

(54) 



with z = Na 2 /7r 2 . The mathematically-inclined reader will immediately notice that the 
coefficients of the first two sums are related to the Catalan numbers. We will postpone till 
later discussion of the evaluation of these sums, in order to focus on obtaining the form of 
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the low-density expansion which is based purely on general considerations. It is clear from 
eq. [50] that in order to have a finite thermodynamic limit, Q{z) must scale as z and 'H(z) 
and I(z) must scale as z 2 for large z. Hence we may define 



111 

lim -Q(z) = g(z) ; lim -^H{z) = f)(z) ; lim -zX(z) = i(z) 

z— >OD Z z— s-oo Z z— >oo Z 



(55) 



where q(z), t)(z) and i(z) have, at most, logarithmic dependence on z. In the limit that N 
and L are large but finite we then have: 



Eo 
N 



2a^p 
M 



(^f) {d - r 2 - Trlog(NA)) 

- 5V 4 + 7rlog(NA)(f) + 2V 2 + Trlog(NA))) + 0(a' 2 : 



.(56) 



This form makes clear that the logarithmic dependence on N must be canceled by q(z), f)(z) 
and i(z) in order to be left with an energy-per-particle that is finite in the thermodynamic 
limit. That is, we have the differential equations, 



d f 

dN 
_d_ 

dN 



(g - V 2 - Trlog(NA)) = ; 



(x + V\ - 5P 4 + 7r log (NA) (t) + 2V 2 + 7r log (NA) 

which are readily integrated to give: 

q(z) = 7T log z + g ; 
f)(z) = —2-k log z + h ; 

i(z) = 7T 2 log 2 Z - 7T (k + 2"P 2 ) log Z + I , 



(57) 



(58) 
(59) 
(60) 



where g, h and i are integration constants. Plugging these functions into eq. [56] we may take 
the thermodynamic limit and we find to 0(a' 2 3 ): 



N 7 



+ (^Y(logV 2 - (21ogAvr 2 + I (2P 2 + &))log 
+ logA7r 2 I(2P 2 + h) + log 2 A7r 2 + £(pf - 5P 4 

+ ^ (« 2 : 



(61) 



There is one further constraint: here we expect that the energy-per-particle should be inde- 
pendent of A up to 0(a' 2 l ) corrections. Using eq. d3 one finds 



A— 
dX { N 



2a' 2 3 p 
Mvr 3 



(2g + h + tt) + 0(a 2 ' 



(62) 
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and therefore there is a further relation between the integration constants, 



2g + h + it = 



(63) 



Note that in eq. EH the energy is completely determined to 0{a' 2 log a' 2 ). Indeed we see 
that the coefficients of the leading logarithms of the form 0(a' 2 n+1 log™ a 2 ) are fixed. The 
change of renormalization scale to obtain a density-dependent coupling introduced terms of 
the form log N, and as the form of the density expansion had to be such as to cancel these 
divergent terms, it is not surprising that the leading logarithms in the expansion can be 
removed by a change of scale. (We will do this explicitly below.) 

In order to go further, one must evaluate the sums, eqs. EH [52] and [53] We evaluate Q 
and H in Appendix II. We recover the form as expected in eq. [60] and find the integration 
constants to be: 



9 
h 



71 ^ 
2 + ^ 
-2tt 



2V- 



2 ■ 



which are of course consistent with eq. [63] We have been unable to determine i 
integration constants g and h then fix the energy density to 0(a' 2 loga^): 



(64) 
(65) 

The two 



£(p) = px 



E 2a' 2 p 2 



N 



M 



1 + (£\ (\oga' 2 - log Xtt 2 + I 



+ 



^ (log 2 4 + 2(l-logA7r 2 )log^ + logA7r 2 (logA7r 2 - 2) + £ (P| - 57> 4 + 



(66) 



It is straightforward to check that this result is in agreement with Refs. [20] and [21] to 
£>(a' 2 3 logc4). As the energy density in the thermodynamic limit cannot depend on the 
geometric constants P 2 and V±, we define 



C 



1 



-An 



5V 4 + i 



While we have been unable to calculate this constant, Ref. [2l| finds 

C A = 2.78 x 10" 3 . 



(67) 



(68) 



In the calculation of Ref. 21 



fluctuations around a mean field 31 



which is based on a systematic EFT computation of quantum 
C arises as a two-loop effect. This is consistent with 
our expectations for the integration constant i, as the leading term in the sum depends on 
the double lattice sums Qo and TZq which are clearly related to two-loop vacuum integrals 
in the continuum limit. 



8 By inspection of eq. I52| it would appear that performing the sum to obtain H. would involve solving 
the two-dimensional double lattice sums Qo and TZ-o in the sense of expressing them as products of one- 
dimensional sums. 
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C. Finite-size corrections 



As we are able to evaluate the sums Q and T-L explicitly, we are able to give the leading and 
next-to- leading finite-size corrections to the thermodynamic limit. Relaxing the thermody- 
namic limit in eq. ESI gives: 



-g(z) = g(z) + ^ + 0(z- 2 ) 

z z 



U(z) 



+ 0{z 



-2\ 



(69) 



In Appendix II we find q = 1 = — f) - Using eq. [50] we see that there are leading and 
subleading 1/N corrections that arise from g , f) ^ 0, and, in addition, there is a correction 
that arises from the binomial coefficient prefactor, as shown explicitly in eq. |50j Taking into 
account both of these contributions gives 



' ^2 



5£ 



FS 



4« 2 P 
MN 



1 + 



(70) 



One easily checks that this expression is independent of A up to 0(a 2 ) corrections. It 
follows that 



5£ 



FS 



£ 



N + °W 



(71) 



which constitutes a ten-percent effect in a system with N = 20 bosons. 



D. Non- universal corrections 

It is straightforward to include the leading non-universal corrections in the energy density. 
By inspection of eq. H7J it is clear that the leading effective range corrections to the ground- 
state energy of N bosons in a finite area vanishes in the thermodynamic limit. Hence, the 
leading non-universal contribution to the energy density is from the three-body force, 



58 



NU 



P 



Do 
36 



(72) 



as one expects on the basis of simple mean field theory considerations. An estimate of the 



leading range corrections to the energy density has been made in Ref. [32 



E. Summary and discussion 

Our final form for the energy density in the thermodynamic limit may be written as 
2a' 2 p 2 



£ + 5£ NU 

2 



1 + (jf) (logc4 - logAvr 2 + I 
+ (^) (log 2 a' 2 + 2(1 - logAvr 2 )log4 + log Avr 2 (log Avr 2 - 2) 



1 



2C 



+ O (a' 2 



+ P 



3 Do 
36 



(73) 
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which is complete to 0(a 2 3 ) if one uses eq. [68]for the constant C from Ref. 2l|, and is valid 



for a 2 loga 2 <C 1. We stress that this result is model independent; the critical assumption 
that we have made is that the potentials via which the bosons interact are of finite range. 
The chemical potential and pressure are readily available from the energy density via 



d£_ 

dp 



V 



-8 + p/i . 



(74) 



The results found in this paper (and also in Ref. H) demonstrate that there is a scale 
ambiguity in the equation-of-state of the Bose gas in two spatial dimensions. While the en- 
ergy density is independent of A and of the renormalization scheme that is used to define the 
coupling constants, in the perturbative expansion this holds only up to 0(a 2 4 ) corrections; 
that is, the inevitable truncation of the perturbative expansion implies that predictions do 
depend on A. In principle, an ideal choice for A will optimize perturbation theory for the 
particular system in question 9 . For instance, if one chooses A = a' 2 /ir 2 , then all logarithms 
of the coupling are absorbed into the coupling itself and one is left with a simple perturbative 
expansion in a' 2 , 



I „2 



£ 



2a 2 P 
M 



71 



I \ 2 
71 



2C 



0(a' 2 



(75) 



where the coupling a 2 is determined self-consistently from eq. [T7] or from eq. [76] below. 

The most natural way of expressing interactions in the EFT is in terms of the Lagrangian 
coefficients, which run with the renormalization group in two spatial dimensions. By con- 
trast, the two-dimensional scattering length is not a natural quantity in the EFT; indeed it 
is the most unnatural quantity that it is possible to form, as it corresponds to the distance 
scale set by the Landau pole. Nevertheless, the energy density can be expressed in terms of 
the two-dimensional scattering length via the formula 



a 2 



71 



log(pA(27r) 2 a|) 



(76) 



which is obtained by comparing eq. [13] and eq. LTHl This is the traditional way of expressing 
the two-body coupling constant [181]. We see that the argument of the logarithm depends 
on A, and is therefore not a physical quantity; any attempt to assign definite meaning to it 
is futile. 

Finally, for facility in comparison, we will express the universal part of the energy density 
in terms of the scattering length. As pointed out in section II, there are various conventions 
used in the literature for the scattering length; one convention, 02, is as given in eq. [18] and 
another identifies the scattering length with the radius of a hard disc, a 10 . In the first 



10 



uiiiesarise in perturbative QCD. Ear a relevant discussion, see Refs. [33( and [3 



Similar scale ami 

Evidently Refs. [32], [35|, [36( claim that Refs. [20( and 2l| are discrepant, and, moreover, that Ref. |2l| is 



incorrect. As pointed out above, we find no discrepancy between these two calculations. We believe that 
confusion may have arisen due to the choice of convention for the two-dimensional scattering length. 
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convention, we have, with A = 1/(2tt) 2 , 



M\\ogpa2\ | log pa^V V 

+ n — - — T lo g 2 I lo g P a l I - 2(1 + log4vr) log | log pa 2 .] 
I log petal V 

+ log 2 4 + log 16(1 + log tt) + log vr(2 + log tt) - 1 - 2C) . (77) 

One readily finds the energy density in the second convention by choosing ci2 = ae 7 /2 and 
A = e^/vr 2 in eq. [76] 

VI. CONCLUSION 

In this paper we have computed the ground-state energy of N identical bosons which interact 
via the most general finite-range potential in a finite area. This energy is expressed as a 
double perturbative expansion in the two-body interaction strength, which is logarithmically 
dependent on the system size L, as well as in inverse powers of L by way of operators that 
break scale invariance at the classical level. Effective range corrections and the leading effect 
of three-body forces enter at 0(L~ 4 ). The structure of the expansion is largely dictated 
by scale invariance and its logarithmic breaking. Indeed, the EFT with the leading two- 
body interaction acts very much like a renormalizable field theory with a coupling constant 
that runs logarithmically. All other interactions beyond the leading two-body interaction 
generate power-law breaking of scale invariance. Using the finite-area ground-state energy 
as a starting point, we have also explicitly evaluated the sums that diverge with powers of 
N and recovered the well-known low density expansion of the ground-state energy density 
in the thermodynamic limit. 

We have seen in this paper that the many-body boson system in two spatial dimensions 
is significantly simpler from a mathematical standpoint than its counterpart in three spatial 
dimensions. The tractability of the two-dimensional system is due both to the logarithmically 
broken scale invariance of the system at leading order in the momentum expansion in the two- 
body sector, as well as due to the expression of two-dimensional lattice sums as products of 
familiar one-dimensional sums. These two features allow one to move smoothly between two 
weakly-coupled quantum regimes that are related by infinite resummations. In particular, 
this tractability allows one to calculate the leading and sub-leading finite-size corrections to 
the thermodynamic limit equation-of-state. In principle, this will enable the quantification of 
finite-size effects in experimental results involving ultra-cold atoms interacting in two spatial 
dimensions. With the results found in this paper, it would be interesting to investigate the 
transition between the confined and thermodynamic-limit regimes using quantum Monte- 
Carlo methods. 

It should be clear that the method presented here for computing the equation-of-state 
and low-density properties of the Bose gas in the thermodynamic limit is not particularly 
efficient. Indeed, the technology developed in Ref. [3lj and carried out in the two-dimensional 
case in Ref. j2lj provides the most efficient and sensible method for treating the low-density 
limit in a model-independent way. Nevertheless, it is interesting to see that the results of 
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the low-density quantum loop expansion can be obtained in an explicit model-independent 
construction without any reference to mean field theory. 



Acknowledgments 

I wish to thank Will Detmold, David Kaplan and Martin Savage for useful discussions and 
for allowing me access to unpublished notes. This work was supported in part by NSF 
CAREER Grant No. PHY-0645570. 



APPENDIX I: Two-dimensional lattice sums 



In an even number of dimensions, it is possible to decompose multidimensional lattice sums 
into products of simple sums (US using methods pioneered by Jacobi |40|. For instance, 
of particular interest to this paper is the sum 



V, 



2.s 



E 

nezVo 



n 



(A-1) 



which is valid for s > 1, where 



CM - E 



m=0 



ft*) = E 



:-i) n 



m=0 



(2m + 1) 



(A-2) 



are the Riemann zeta function and Dirichlet beta function, respectively. The case s — 1 
that we are interested in is singular as ((s) has a simple pole at s — 1. In order to subtract 
this pole, we require the Laurent expansions of £(s) and (3(s) about s = 1 [4j]]. We have 



COO 

0(8) 



s-1 

71 



+ 7 + 0{s-l) ; 



l + ( S -l)( 7 + 21og2 + 31og7r-41o g r(i))+ O ((s - l) 2 ) , (A-3) 



where 7 is Euler's constant. One then easily finds 

V-2 = lini 1 > — — — j = 47rlog(e2 7r-3r(|) 



S-1 



(A-4) 



It is straightforward to generalize eq JA-ll to 

1 A (-1! 



E 



(n 2 - 17 2 ; 



*E 



(2£+V 



C[s,l 



rj 



(2£+l) 



(A-5) 



which is valid for s > 1. Here ((s,a) is the Hurwitz zeta function [41]. Using methods 
similar to those outlined above, one readily finds 



5} (E 

Wo ^" 



7T 



S-1 



P 2 -7T 7 -4^ 



£=0 



(-1)* 

(2£ + r 



(2£ + i; 



, (A-6) 
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where t/>o is the digamma function. Finally, one can write 

M ,) = -_ +?! - I7 _4g±J_ ^ (i _ _R_ j . (A . 7) 

APPENDIX II: Catalan sums 

In this appendix we evaluate the sums which diverge with powers of N in the thermodynamic 
limit. The first sum we wish to evaluate, eq. EH may be expressed as 

OO 

= ^(-^) n C(n)z n V 2n , (A-8) 

n=2 

where the C(n) are the Catalan numbers 11 , which have the integral representation 42| 

4 f°° uj 2 ( 4 \ n 

7T J (1 + W 2 ) 2 Vl 1 +^ )/ 

Using eq. IA-11 we can write 

ew = t T l££r'E g ' ( < nWn) • (A ' 10) 

70 ^ ' n=2 

where z = — 4z/ (1 + w 2 ). By expanding eq. IA-61 and comparing with eq. IA-11 it is straight- 
forward to find 

t rmm - - i&rM'-w^))- (A " U) 



Using the asymptotic form of the digamma function for large argument as well as the Dirich- 



let sums 41 



one finds 



G(~) = ;r: (logz+^ + ^) + 1 + 0{z~ x ) . (A- 13) 



And finally, matching to eq. [691 



q(z) = vrlogz + | + V 2 ; fl = 1 • (A-14) 



11 C(n — 2) is the number of ways in which a regular n-gon be divided into n — 2 triangles if different 
orientations are counted separately [41] (Euler's polygon division problem). They are related to the 

2n\ 



central binomial coefficients via C(n) = ( ™)/{n + 1) 
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Similarly, eq. [53] may be written as 

oo 

U (*) = E(- 1 ) n+1 (" + VC(n)z n+1 V 2n . (A-15) 

n=2 

Using the integral representation [4l[ 



(B+1)C(B) = ^Cjrh)(iT^^ d "' (A - 16) 



7T 



and proceeding as above one finds 

H{z) = -2nz 2 Mogz + l + ^J - z + O(z ) . (A-17) 

And finally, 

= -2vrlog2 - 2vr - 2V 2 ; f)o = -1 • (A-18) 

As the convergence properties of the Catalan sums determine the transition to the low 
density regime, we can get a more accurate measure of the expected region of validity of 
eq. 1471 through a simple convergence test. According to DAlembert's ratio test, the sum, 
eq. \A-8\ is convergent when 

lim — ni Yd z < 1 ' ( A ~ 19 

from which one easily finds z < 1/4, or Na^ < 7r 2 /4. The same estimate follows from the 
other sum, eq. IA-151 
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